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Abstract — The paper describes an impedance measurement 
technique using a particular sampling method which is a time 
saving and helpful method as compared to the complexity of to 
DFT calculation. In this technique a sinusoidal signal is used as 
an excitation signal and the response of the signal that are 
proportional to current flowing through and voltage across the 
measured impedance got after by sampling are extracted and 
processed. The impedance of the object under consideration is 
calculated without using Fourier transform. This method will be 
first evaluated in MATLAB by means of simulation. 

This paper is definitely a comparison of two impedance 
measurement methods namely the DFT and particular sampling 
in the domains of error, time savings and standard deviation. 
The Labview application is used to generate the excitation signal 
using the data acquisition card NI-USB 6251. The sampled 
responses are compared by the graphical program developed in 
the labview platform is used to compute the error, standard 
deviation and percentage of time saving. 

Index Terms — Data acquisitions, Fourier transform 
Impedance, Matlab, Labview 

LIST OF IMPORTANT SYMBOLS ABBREVIATIONS 


AC 

- Alternating Current 

dB 

- Decibels 

DC 

- Direct Current 

DDS - 

Direct Digital Synthesis 

DFT 

- Discrete Fourier Transformation 

DSP 

- Digital Signal Processing 

e 

- Electrons 

EDF 

- Electricite de France 

F 

- Farad 

f 

- Frequency 

FFT 

- Fast Fourier transformation 

G 

- Conductance 

Hz 

- Hertz (unit of frequency) 

I 

- Current 

IET 

- Institution of Engineering 

IEEE 

- Institute of Electrical and Electronic Engineers 

k 

- Kilo 

LCR 

- Inductance (L), Capacitance (C) and Resistance 

M 

- Mega 

n 

- Nano 

P 

- Pico 

PSoC 

- Programmable System on Chips 

S 

- Siemens(mutual inductance) 

SIR 

- Source Impedance Ratio 

SST 

- Single sine technique 

XC 

- Reactance of Capacitance 

XL 

- Reactance of Inductance 

R 

- Resistance 

Z 

- Impedance 

Q 

- Ohms (Unit of reactance) 

D - 

Angular frequency 

9 

- Phase angle 
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I. Introduction 

The terms impedance and resistance both means oppose to the 
flow of current. In circuits working on direct current, only 
resistors produce this effect. On the other hand in alternating 
current (ac) circuits, the other components, like inductors and 
capacitors, also provides opposition to the flow of current. 
The total impedance in a circuit is the sum of opposition of all 
the elements to the flow of current. The opposition provided 
by inductors and capacitors in a circuit is called by the same 
name reactance, characterized by X and calculated in ohms 
(Q). Since the symbol for capacitance is C, capacitive 
reactance is represented by XC. Similarly, since the symbol 
for inductance is L, reactance offered by inductance is 
represented by XL. Inductors and capacitors not only affect 
the magnitude of an alternating current but also its time 
dependent properties - or phase. When most of reactance to 
current flow in a circuit is offered by a capacitor , a circuit is 
said to be largely capacitive and the current through the circuit 
is said to be leading the voltage in phase angle. When most of 
the reactance to the current flowing through the circuit comes 
from an inductor, a circuit is said to be largely introductory 
and the current is said to be lagging behind, the applied 
voltage in phase angle. The more introductory a circuit is, the 
phase difference between the current and voltage approaches 
90 degrees. It’s sometimes easier to perform calculations 
using admittance, the mutual impedance. Admittance is 
characterized by Y and calculated in Siemens (S). Like 
impedance, admittance can be explicate as a intricate figure, 
where the conductance, the mutual of resistance, is the real 
item, and the susceptance, the mutual of reactance, is the 
imaginary item. 

There are numerous specialized and organic objects, the 
parameters of which can be evaluated by impedance 
estimation. A case of such approach is checking and 
diagnostics of anticorrosion protection of huge specialized 
items like pipelines, extensions, fuel tanks and so on. A case 
of natural security and wellbeing is the utilization of 
impedance estimation for observing of water permeation of 
dams. Another illustration is detecting innovation e.g. relative 
humidity sensors. In the event of the initial two cases there is a 
need of circulated estimations in numerous spots - various 
sensors are required, so the estimation time is expanded 
definitively or the required number of instruments working in 
parallel gets to be higher. If there should be an occurrence of 
sensors, the aggregate cost of the gadget, circuit 
measurements and power utilization ought to be kept as little 
as would be prudent for the most part because of monetary 
reasons. The above displayed rules prompt scanning for new 
strategies which permit to disentangle impedance estimation 
instrument. 
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Local impedance calculation methods are always revised 
because of highly effective digital and analog converters 
development, making use of the digital signal processing 
methods and fast construction of recent and more energized 
digital executing signals which can effectively undergo a lot 
of operations in each second. Measurement of impedance is a 
very essential subject with various recent findings researched 
newly. Few research works have been targeted on 
measurement of the impedance sensors while some utilize the 
newly transmitted network executing algorithms to enhance 
impedance calculation. An easy impedance calculation 
technique staged and improved minimized the front-end 
analog circuit to diminish the effect of 
frequency dependence. This method is staged with 
volt-ampere method with anonymous impedance and 
reference impedance arranged in series and delivered with 
same electric current from a sine generator. Voltage across 
each impedance is acquired simultaneously with dual analog 
to digital channel converter with the phase of impedance 
extracted with aid of signal processing digital algorithm. 
Algorithms with sine-fit structure have an option to calculate 
the signal parameters [15] within a group of acquired samples. 
Ever since gap within the sampling rate and sine wave 
frequency mostly are not known by accuracy, the said 
algorithm should be able to calculate sine wave frequency as 
well. It allows the issue of regression of nonlinear which is 
sorted by using procedure of iterative that is named sine-fit 
four parameters algorithm. Take for instance the 
measurements for impedance, the structure algorithm may be 
used on the channels individually, moreover the final 
researched probability may be diminished by moving the 
common rate of occurrence with force in a seven sine-fit 
parameter algorithm [15]. Nevertheless, the iterative 
impedance structure of the sine fit output large amount of 
processes that must be executed. Another problem that may 
hinder effectiveness of the algorithm developed with device 
such as DSP is memory management and specifications. It 
should be noted that the seven sine-fit algorithm developed 
2N rows and seven column matrix, in accordance to the sum 
of acquired samples by both channels. 

This large matrix and its manipulation together with the 
limited memory available in the DSP restrict the number of 
samples that can be processed. In this paper, the algorithm 
proposed is shown to be applicable in a memory restricted 
system by bypassing the need to build the 2N*7 matrix. This 
modification allows DSP to process more samples to calculate 
the sine wave frequency that must be accurately known for a 
correct calculation of the sine signals’ amplitudes, phases and 
DC items by the sine-fit algorithm [15], in impedance 
measurements an accurate frequency value is not required. 
The time information, and thus the sine signals frequency can 
be removed from the problem by making an XY plot of the 
two waves. The result is an ellipse whose values can be 
calculated by ellipse-fit algorithm such as the one presented 
and then improved. The sine signals values can then be 
extracted from the ellipse parameter as proposed. This 
algorithm has been recently optimized and adapted for use in 
DSP based impedance measurements. In the modified 
ellipse-fit method published [15], there is no limit on the 
number of samples that can be used since only nine values 
must be stored and the samples themselves can be discarded 
after their contribution to the nine stored values is taken into 
account. These properties make the algorithm a prime 


candidate for efficient in a DSP based impedance 
measurement instrument. 

The aim of this master diploma thesis is to present and 
implement impedance measurement method using particular 
sampling method [13], which is a substitute for DFT 
calculation. The method uses sine excitation signal and 
sampling response signal proportional to current and voltage, 
flowing through and across the measured impedance. Fourier 
transform is not used in this method. The method will be first 
assessed in MATLAB by means of simulation. The excitation 
signal is generated using National Instruments data card and 
graphical software platform LAB VIEW. The software is 
developed to measure the voltage across the measured 
impedance and current flowing through it and hence calculate 
the impedance and verify it with the conventional DFT 
method. 

The second chapter describes various methods of impedance 
measurements such as null method using different kinds of 
bench and also different meters workings on the basis of 
resonance method and active method. The third chapter deals 
with the theoretical aspects like phase-sensitive detection, 
DFT, FFT, sine-fit and ellipse-fit algorithm and Goertzel’s 
algorithm, and detailed explanation of complexity of each 
method.. 

The next chapter explains the concept of particular sampling 
and also explains the result of simulation using MATLAB and 
the implementation of software part for the impedance 
measurement using LAB VIEW for the particular sampling 
method. This chapter deals with the experimental data 
analysis on a circuit and trying to find the impedance using 
DFT as well as particular sampling method and determining 
the error in impedance measurement The next chapter 
describes the experimental results of the software based 
implementation of the particular sampling method. The next 
chapter includes the results based on the experimental results 
and scope of improving the results and thus reducing the 
complexity, the chapter also deals with, what have been done 
in the field of diploma thesis. The document is finished with 
bibliography section and list of figures and tables. 


II. IMPEDANCE MEASUREMENT METHODS 


The common methods used for impedance measurements 
are: auto balancing bridge method, resonant method, network 
analysis method. [20]. 

The impedance is used to specialize electronic circuits and 
items. At given frequency when it is passed through electronic 
device or circuit impedance is defined as opposition to AC 
(alternating current) 

+j 

1 Z = R+jX = |Z|,:0 


Z{R, X) 



Real axis 


R = jZ| cost 1 
X = |Z| sinfl 

f |Z|= * +x 2 

1 0 = tali' (*) 


Fig. 2.1. Impedance vector analysis [10]. 


Fig. 2.1 shows that the impedance has an imaginary reactance 
X and real part resistance R. In order to measure impedance, 
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two values are required to be calculated as it is a complex 
quantity. Both real part and imaginary part of impedance 
vector are usually used while measuring all the impedance 
measurement instruments. Parameters such as |Z|, 0, |Y|, R, X, 
G, B, C and L, later on instrument converts these real and 
imaginary parts. 

There are various impedance measurement methods based 
on various parameters, such as measurement range, frequency 
of operation, ease of use and measurement accuracy. The 
impedance measurement methods like resonant method, 
bridge method, I-V method, auto balancing bridge method 
and active methods are explained in the following sections 
with their advantages and disadvantages. As the recent 
developments in impedance measurement devices LCR 
meters are more commonly used for the accurate computation 
of impedance 

In order to measure impedance, the single sine technique 
(SST) is commonly used. By repeating the measurements at 
various frequencies, the impedance spectrum can be obtained 
directly from measurement results as a function of frequency 
in a range of a few decades. SST is based on excitation of the 
object with a harmonic signal and vector measurement of two 
signals: voltage across and current through the tested object. 

In the digital implementation of the impedance 
measurement method excitation signal is produced with the 
aid of Direct Digital Synthesis (DDS) using a D/A converter 
and memory containing sine samples. To extract signals 
proportional to the voltage across (Vy — Vx ) and current 
through (lx — ly ) the calculated impedance Z the input 
circuitry has to be used. The block diagram of an impedance 
measurement system is shown in Fig. 2.2. The construction of 
the input circuit is very important as parasitic capacitances 
and real-life parameters of the operational amplifiers can 
significantly influence a measurement result. Signals and are 
sampled synchronously with the clock generator using two 
A/D converters and placed in memory in the form of two sets 
of N samples of signals, Vy [n] and Vy [n] 



Fig. 2.2. Block diagram of the impedance measurement 
system. 


To determine the orthogonal parts of the acquired signals 
on the basis of the collected samples, various techniques can 
be used: the use of Fourier transformation (DFT or FFT), 
sine-fitting algorithms, ellipse-fitting algorithms and others. 

The methods of computing impedance describe the 
standard ways of computing impedance in the electric circuit 
using the concepts such as 


•Null method 
•Resonant method 
•Active method 

The strength and weakness of each measurement method 
depend on accuracy, cost and effectiveness. Sometimes more 
than one method is used to compute the impedance in the 
electric circuit which is dependent on the characteristics of the 
circuit 

2.1. Null method 

The principle of null method is the Wheatstone’s bridge 
and all instruments will have an appearance as the same. The 
voltage is applied on one of the portion of circuit and the 
parameters of the circuit are adjusted to obtain the balancing 
condition and at this moment there exist a mathematical 
relation between the known and unknown components in the 
circuit and by which the unknown impedance can be 
calculated 



Fig. 2.3. Wheatstone’s bridge [1] 


Fig. 2.3. represents a Wheatstone’s bridge. In the circuit 
diagram Rx is connected to the unknown terminals and the 
resistance Rl, R2, R3 is connected to the known terminals. 
The resistance connected in the known terminals can be 
adjusted to achieve the balanced condition so the 
galvanometer (G) connected shows the null deflection and 
thus the nodes B and C are at same potential. 

In the balanced condition we have the following 
mathematical equation [1] 


R3 Rl 
Rx _ R2 


So we have 


R2 

Rx = — * R3 
Rl 


( 2 . 1 ) 

( 2 . 2 ) 


Thus the unknown impedance can be calculated from the 
above equation (2.2). For the ac circuits the resistances are 
replaced by admittances or impedance, battery by a signal 
generator and the galvanometer by an ac voltage detector. 

The major advantages of the null method are the accuracy 
of the measurement. The accuracy is dependent on the factors 
like dial resolution, internal shielding and constructional 
factors. Most null detectors are dependent on signal frequency 
and are designed to operate with high accuracy at an ambient 
frequency, below this frequency the accuracy is decreased due 
to the sensitivity and above this frequency the accuracy is 
diminished due to the residual impedance. 

The different kinds of null instruments are illustrated below 


> Maxwell’s bridge 
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B 



Fig. 2.4. Maxwell’s bridge [1] 

Fig. 2.4. shows the Maxwell’s bridge, the balancing 
equations [1] are as follows, 


Rx = R2 * R3 * AGs 

(2.3) 

Lx = R2 « R3 * AGs 
Where Rx and Lx 

(2.4) 

are the unknown resistance and 


inductance. This bridge measures the inductance and 
resistance connected series in terms of standard capacitance 
(Cs) and conductance standard (Gs). This bridge is ideal for 
measuring low impedance at high frequencies. The initial 
balance is achieved respect to the short circuit placed at the 
unknown terminal. 

Advantages 

•High accuracy 

•Best among other instruments to measure very low 
values of impedance and series resistance. 

•Broad measuring range of resistance 
Disadvantages 

•Maximum frequency limit is 30MHz 

•Difficult to calibrate 

•Slow 

> Schering bridge 

The Schering bridge is also called radio arm bridge, and is 
shown in Fig. 2.5. It is specialized to perform more specific 
measurements in one form or another and is widely used in RF 
bridge circuits. 


B 



The balancing equations [1] are as follows 


R1 R2 R4 



(2.5) 

C4 “ C3 “ Cl 

R2 


Rx = — *AC1 

C3 

(2.6) 

Cx = AC4 

(2.7) 


Rx and Cx are the unknown resistance and capacitance. 
This bridge measure admittance in terms of parallel resistance 
equivalent and either positive or negative equivalent parallel 
capacitance or thus both capacitive and inductive impedance 
can be measured. Initially the balancing condition is achieved 
with respect to an open circuit and first order errors due to 
signal generator are not present. 

Advantages 

•Broad frequency range 

•Less calculations needed 

•Broad range of resistance value measurement 

•Good accuracy 

Disadvantages 

•Minimum value of frequency is 500kHz 
•Operator training required 

•Small range of inductance and capacitance 

measurement 

•Slow 

> Admittance ratio bridge 




Fig. 2.6. Admittance ratio bridge [1]. 

Fig. 2.6. represents an admittance ratio bridge, the balance 
condition equations [1] is as follows, 


R1 

(2.S) 

Gx = — * Gs 

R2 

R! 


Cx = — * Cs 

R2 

(2.9) 


Gx and Cx are the unknown conductance and capacitance. 
This bridge has its application in the low frequencies and 
extends up to 1 MHz. The circuit can be used as a comparator 
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when external standards are available. Initial balance is made 
with respect to an open-circuit. 

Advantages 

•Wide range of capacitances 
•Good accuracy 

•Can be used as comparator using parallel range 
extension 

Disadvantages 

•Cannot be used to measure small impedances 
•Limited to frequency below 5 MHz 
•Operator training needed 

> Transformer ratio bridge 




Fig. 2.7. Transformer bridge [1] 


Fig. 2.7. presents the transformer bridge and the balancing 
condition equations [1] are, 



/Nx 

nx\ 


(2.10) 

Gx = 

( 77 “ =■ 

* — ) 

* Gs 


\Ns 

ns/ 



/N'x 

nx\ 



Cx = 

1n7 

H 

ns/ 

* Cs 

(2.11) 


The Gx and Cx are the unknown conductance and 
capacitance. This is a very special bridge that offers very wide 
range of measurement values than that offered by 
conventional Wheatstone bridge based bridges. Primary and 
secondary turns-ratio in in both the input and te detector 
positions are tapped to provide the effect of variable internal 
standard, thus eliminating the need of adjustable resistors and 
capacitors , but they can be added for higher resolutions. 
Advantages 

•Wide range of frequency 

•Wide measurement range for resistance and 
inductance or capacitors 
•Can be automated 

•Can be made grounded, balanced or unbalanced 
based on measurement requirements 
Disadvantages 

•Very slow 

•Requires skilled operator. 

> Twin -T bridge 




Fig. 2.8. Twin-T bridge [1] 


Fig. 2.8. presents a Twin-T bridge and the balancing 
equations [1] are 



Cx = AC2 (2.13) 

The Gx and Cx are the unknown admittance and 
capacitance. The twin-t bridge is actually an admittance 
bridge which offers a very high accuracy. The initial balance 
is made by means of an open circuit and the main 
characteristics of the bridge is due the fact that bridge balance 
is frequency sensitive, represented by the factor oj 2 in the 
balancing equation (2.13). 

Advantages 

•Wide range of frequencies 

•Wide range of impedances of capacitance, 
inductance and conductance 

•Can measure conductance in terms of capacitance 
increment 

•Most accurate method for measuring the 
conductance 
Disadvantages 

•Expensive 

•Slow 

•Cannot measure small resistance 
•Requires skilled operator. 

> Q bridge 


B 
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Fig. 2.9. presents the Q-bridge and the bridge equations [1] 
are as follows 

L1 , , 
Lx = — * ACS (2,143 



The Q-bridge actually works on two principles, a situation 
made possible through the use of a phase sensitive detector. 
The unknown and known components are arranged in the 
Wheatstone’s configuration, and the detector separated the 
unbalanced voltage between the points A and C into the 
in-phase and quadrature elements. The in-phase elements are 
due to the reactance unbalance. The balancing between the 
reactance is achieved by adjusting the in-phase unbalance 
voltage to the generator reference voltage 8,;. 

Advantages 

•Has very range of measuring 
•Versatile- can be used to for a wide variety of 
components 
Disadvantages 

•Low accuracy 
•Slow 

•Measuring range of Q is from 5 to 500 
> Thurston bridge 


SLIDING 

SHORT 



Fig. 2.10. Thurston bridge [1]. 

The Thurston bridge is presented in Fig. 2.10 and is a type 
of microwave bridge, it can be used for wide range of 
frequencies. The voltages are induced into the junction of the 
bridge branches, by means of conductive loops that are 
inductive in nature and the voltage induced is proportional to 
the angle of orientation of the plane of loops, with the center 
of conductors in the respect branch [1]. 

I = kVY (2,16) 

In equation (2.16), V is the generator voltage Fig. 2.10, Y is 
the branch admittance and k is the coupling coefficient, which 
is proportional to the angle of orientation of plane of the loop 
with center of the conductors. 

Advantages 

•Suitable only for measurements above 40MHz. 

•Can be used to measure transistor parameters. 

•Commercially available. 


Disadvantages 

•Low accuracy and resolution. 

•Very slow. 

•Cannot be used for measurements at low frequency. 
> Young bridge 



Fig. 2.11. presents the Young bridge and the balance 


condition equations [1] are 


G x = kl. G^ — G^J 

(2.17) 

Q: = k(C A — C c J 

(2. IS} 


Where in the equations (2.17) and (2.18) G x , are the 
unknown conductance and capacitance. 

This device is specially designed to measure small values 
of the capacitance. The working is based on the conversion 
principle of wye configuration [1] of admittance between the 
corners of the bridge A, B, in Fig. 2.11., and the capacitance is 
equivalent to the differential capacitor and conductance 
standards to maintain the constant and equal admittance in 
arms AB and BC Fig. 2.11. 

Advantages 

• Very useful for measuring very small 

admittances. 

• Good accuracy. 

Disadvantages 

• Can be used for used only at fixed frequencies 

• Higher frequency limit is 1 MHz 

• Is slow. 

2.2. Resonance method 

Resonance methods are used when the Q value that is the 
ratio of reactance to the resistance is very high and are either 
connected series or parallel. Depending upon the arrangement 
either an ammeter or a voltmeter will be used as a detector. In 
both case unknown component value is determined by 
changing the frequency, or resistance or the reactance. In most 
case the variation of frequency is not a good option for high 
accuracy measurements because the frequency response 
versus the voltage (or current) is not symmetrical about the 
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resonant frequency so it can lead to errors. This method is 
ideal for three terminal measurements and cannot be used at 
frequencies that are closer to self-resonant frequency of any 
individual component. In series resonance measurement, the 
frequency and voltage are held constant. The circuit diagram 
is shown below 


a) 




RESISTANCE, REACTANCE 
OR FREQUENCY 


Fig. 2.12. Series resonance [1] 
a) Circuit, b) Frequency response 
The series resonance circuit and frequency response is 
shown in Fig. 2.12. 

In all cases qj L = — and for resistor variation we have the 

■iiiC 

following result [1], 


Rs = 


^-i 


-* R 


( 2 - 19 ) 


In 


The standard resistance, R, is removed and added in the 
circuit to lower the values of the resistance, Rs of the series 
circuit. 

For reactance variations, 


Rs = 


utL-L,) 



For frequency variations, 


Rs = 


LU) - LOjJ 



( 2 . 20 ) 


( 2 . 21 ) 


In parallel resonance measurements, variations in 
conductance, susceptance or frequency variation is used, 
a) b) 




SUSCEPTANCE 
OR FREQUENCY 


Fig. 2.13. Parallel resonance [1] 
a) Circuit, b) Frequency response 


observed as a means of standard conductance G, in and out of 
the circuit. 

For conductance variation 


Gp = 



For variation in susceptance 


Gp = 


UC - O 



For frequency variations, 


Gp = 


rttO — £i5 r J 



( 2 . 22 ) 


(2.23) 


(2.24) 


The different kinds of resonance instruments are described 
below 

> Q- Meter 



Fig. 2.14 Q-meter [1] 


The Q-meter is presented in the Fig. 2.14 and is the widely 
used impedance measurement unit among the other available 
apparatus. It is available in many forms, the most common 
version reads directly in Q and the value of resonating 
capacitance, at particular frequencies in a direct visible scale 
of measurement. Q-meters which are equipped with a 
resistance used for the insertion of the test signal voltage 
cannot be used for making measurements across the 
components with external dc bias applied because of the risk 
of damage, due to external resistor or the thermo couple unit, 
or even the bias can affect the calibration of Q. 

Measurement equations [1] 



Vr 

VI 


(2.25) 


Fig. 2.13.a) presents the parallel resonating circuit and b) 
represents the frequency response of the circuit. In all cases, 
(jjL = —, In this method, the frequency and current are kept 

CiL 

constant during the measurement. Voltages V rl and V r2 are 


4: = ~—p - (2.26) 

* C r 

Advantages 
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•Easy to carry out measurements at different 
frequencies. 

•Fast operation as compared to manual bridges. 
•Wide frequency range. 

Disadvantages 

•Frequency monitoring is needed for accuracy. 
•Uncertainty in Q due to high frequencies. 

> Immittance trans-comparator 


— * 

_ V 

Z=| (2.29) 

From the definition, active methods of measurement are 
defined as the measure of the ratio of complex voltage to 
complex current following the above condition stated by 
Ohm’s law. 



b) 



FREQUENCY OR 
CAPACITANCE 


Fig. 2.15. Immittance trans-comparator [1] 
a) Circuit diagram b) Frequency response 


Fig. 2.15.a) presents an immittance trans-comparator 
circuit and b) represents the frequency response of the circuit. 
This equipment works on the principle of parallel resonance 
or anti-resonance. This can be used even for the measurement 
of di-electric constant. The machine build in various forms 
can be used either by changing the frequency or capacitance 
to measure the Q or the bandwidth of the parallel resonant 
circuit. 

Measurement equations [1] are 


G = 


AC 


f ± 

Af 


(2.27) 


1 

mL = —— (2.20) 

ojC 

This is not a field instrument as it requires a highly stable 
and tunable source with more power capability. 

Advantages 

•Accuracy. 

•Can be used for a wide range of frequencies and 
even can be used to measure positive and negative 
impedance. 

Disadvantages 

•Very slow operation. 

•Need skilled operator. 

•Capable only for measuring high Q value. 


2.3. Active method 

The actual definition of impedance is by Ohm’s law, is the 
ratio of the complex voltage to the complex current 


* , , IV| < & v |VI 

Z = |Z| < 0 = |r| ; 0 = -^ < 0 v -< Q (2.30) 

Where I Z|, |V|, HI are the absolute magnitude values of 
the impedance, voltage and the current and the phase angles 
are represented by 0 V , 0j respectively. The different 
kinds of active instruments are: 

> Vector impedance meter 



ALC 

Fig. 2.16. Vector impedance meter [1] 


The vector impedance meter is shown in Fig. 2.16. 

Measurement equations is as follows 

Z = f (2.31) 

These are the most recent devices and are convenient over 
the conventional bridges. Co-axial probes are used to connect 
the unknown component whose impedance is to be measured 
to the measuring circuit, and impedance magnitude and the 
phase angle can be read on the digital meter. The major 
applications are in testing, laboratories and also for design 
purposes. 


Advantages 

•Faster than bridges. 

•No prior skills needed to operate. 

•Versatile in nature, can be used for a wide range of 
frequencies. 

•Wide range of measurements. 

Disadvantages 

•Low accuracy. 

•Possibility of error at very high frequencies. 


> L-C meter 


18 


www.erpublication.org 
























International Journal of Engineering and Technical Research (IJETR) 
ISSN: 2321-0869 (O) 2454-4698 (P) Volume-8, Issue-4, April 2018 



Fig. 2.17. L-C meter [1] 


An LC-meter is presented in Fig. 2.17. Measurement 
equations [1] are as follows. 


Vo 

(2.32) 


Vo 

^T«C (2.33) 


input signal V s (t) 



reference signal V r (t) 


lock-in 

amplifier 


amplitude 

phase 


Fig. 3.1. Lock-in amplifier [2] 



DUT ■ 

w 




sine wave 
generator 


V s (t) 


V r (t) 


+ lock-in 
> amplifier 


R 


0 


The L-C meter can only be used for measurement at fixed 
frequencies and does not provide any information about the 
resistance or conductance of tested reactance. The L-C meter 
sometimes is equipped to a dc output which is proportional to 
the inductance and capacitance that is measured and can be 
used to plot an XY graph using a plotter. 

Advantages 

•Faster measurement 
•No operator experience needed 
•Good accuracy 
•Portability 

•Linear dc output for data processing 
•Easy to calibrate and is inexpensive. 

Disadvantages 

•Operates only in fixed frequencies 

•Limited capacitance and inductance measurement. 

III. IMPEDANCE ORTHOGONAL PART 
DETERMINATION METHODS 

In this section, the different theoretical concepts and 
algorithms for impedance measurements are explained. The 
different DSP concepts like phase-sensitive detection, DFT, 
FFT, sine-fit algorithm, ellipse-fit algorithm, Goertzel filter 
and particular sampling are explained 

3.1. Phase-sensitive detection 

Phase-sensitive detection as shown in Fig. 3.1. and Fig. 3.2. 
is used to acquire very small signals in the presence of large 
additive noise. This concept led to the development of lock-in 
amplifiers, which makes this technique possible to execute. 
Lock-in amplifiers use the fact about time dependencies of a 
signal to extract it from a noise accumulated signal. A lock-in 
amplifier performs a signal multiplication to the input with a 
reference signal, this process is also called heterodyne 
detection, and then applies a low-pass filter to the result as 
shown in Fig. 3.3. This is called demodulation or 
phase-sensitive detection and recovers the signal the 
particular frequency that is in area of interest [3]. 


Fig. 3.2. Concept of lock-in amplifier [2] 


mixer LP filter 



mixer LP filter 


Fig. 3.3. Phase sensitive detection [2] 

The reference signal is either generated by the lock-in 
amplifier or provided externally from a signal generator. The 
phase sensitive detection can be generally described as 6 stage 
process and they are described as follows [17]. 

• Modulation - The input signal is modulated or mixed up to 
a certain frequency 

• Pre-amplification - A high-speed amplifier is used to 
amplify the signal to a desired level to feed to the circuit. 
Sometimes the amplification is also proceeded by a high pass 
filter in-order to remove the DC components in the signal. 

• Reference signal - In most cases a pulse signal at the 
modulating frequency is used, it can be from a function 
generator or astable multi-vibrator. 

• Multiplier - In this stage the pre-amplified signal is 
multiplied with the reference signal. The more generally used 
multiplier is an analog switch. The reference signal is used to 
open the analog switch periodically and letting the modulated 
input signal to pass through it. 

• Integrator - The multiplied signal is passed through an 
integrator. Since the multiplied signals have a lot of frequency 
components in it, all of them will change to zero after 
integration, but the component that have the same phase and 
frequency remains non zero, and this component corresponds 
to the product of modulated signal with reference signal. 
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• Low pass filter - After the integration stage, the signal is 
demodulated to DC output. Due to presence of high frequency 
components we use low pass filter to filter the signal at certain 
frequency range. 

Suppose fl and f2 be two signal characterized by the 
following wave equation fl = A sin to L t 

andf2 = Bsinto 2 t. 

Now f 1 *f2 is given by 

AB sintL^t sin(td 2 t T tp) = — — * {cos[( -f to 2 ) * 
t -|- cp] — cos[(to L — to 2 ) * t — tp]} (3.1) 

Let’s integrate the product over the time [0, T] 

if T ir T ab 

- 1 fl *f2dt = - 1 -— * {cos[( ^ + q 2 ) * t+ <p] - m[( u 1 - (i3 2 ) * t - tp]} dt (3.2) 

T Jfl 1 Jn 2 



Fig. 3.5. Computation of real and imaginary part using correlation 


The mathematical equation to compute DFT of a signal is 
described below, 

For a discrete time signal f(n), the N point DFT is defined as 


AB [ T AB f T 

= -— J ™[(M“i)*t+<p]dt+—j ( 3 . 3 ) 

When T botl>t he integral parts will be zero, except when 
_ AB 

iii L = and sp = 0 the integration result is given by— 

3.2. Discrete Fourier transform and Fast Fourier transform 


—i -iarmk 

F[k] = ) f(n)e H ,k = 0.1.2. . . . . . N - 1 ( 3.4) 

H-U 

The matrix method to determine the DFT of the discrete 
time signal is given below 


DFT is the mathematical process of extracting the 
frequency components from a time domain signal. The DFT 
operates by comparing or correlating the signal to be analyzed 
against sinusoidal waveform or sinusoidal basis functions. 
The comparison is achieved by the mathematical operation 
called correlation as depicted in Fig. 3.4. 

Correlation is either a measure of similarity of the signals 
or how strong the signal is present in the other. The analysis 
base signals used should be of the same length of the signal 
under processing. The sinusoidal basis functions are both sine 
waveforms and cosine waveforms. The comparison of the 
signal being analyzed with the sinusoidal basis functions 
results in complex numbers. The result of correlation of the 
signal with the sine wave sinusoidal 

basis function is stored as the imaginary value and the result 
of correlation of the signal being analyzed with the cosine 
wave sinusoidal basis function is stored as the real part as 
shown in Fig. 3.5. 

The magnitude of the results of correlation is used to 
compute the magnitude spectrum. And from the complex 
number result the phase values are computed to obtain the 
phase spectrum. 


Fig. 3.4. Schematic showing working of DFT 
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functions 


Signal being analysed 
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(3.5) 


(3.6) 


In the equation (3.6), f[0], f[l], f[2].f[N-l] represents the 

discrete time signal and F[0],F[1],.F[N-1] represents the 

corresponding DFT’s and in the matrix W is called the 
twiddle factor matrix or matrix of linear transformation and 
given by 


W = e 



and W = W !H 


(3.7) 


The FFT is a fast algorithm used for computing the DFT. In 
FFT algorithm we divide calculations to compute 2-point, 
4-point DFT’s and generalize to 8 -point, 

16-point,., 2 r p oint. 

To compute the DFT of an N-point sequence using DFT 
equation (3.4) would take ft 11 multiplications and additions. 
The FFT algorithm computes the DFT using Nlog z N 

multiplies and adds. 

When we consider the cost of computation FFT is faster 
than DFT for big values of N To compute the DFT of an 
N-point sequence using the equation (3.4), would require N z 
complex multiplies and adds, which is together 4N 11 complex 

H 

multiplications and additions. There are 7 butterflies per 

stage, and \og 2 N stages, so that means about 
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4 * - * log 2 N = 2E^ log 2 ft real multiplies and 3Nlog s N 
real adds for an N-point FFT. 


x (i) = [[D®] T D (il ] _ 1 [[D t 0 ] r y] 


(3.9) 


Table 3.1. Cost of calculation comparison of DFT and FFT 


N 

Total DFT 
calculations 

4N 2 

Total FFT 
calculations 
2Nlog z N 

Speedup 

2 

16 

4 

4 

4 

64 

16 

4 

8 

256 

48 

5 

1024 

419304 

20480 

205 

65536 

1.7 * 10 1C 

2.1 * 10 6 

10 4 


3.3. Sine-fit algorithm 

Sine-fit algorithms are an option to estimate the sine-signal 
parameters from a set of acquired samples. Sine-fit algorithms 
were standardized for the characterization of ADCs [11] [12]. 
There are three parameter-sine fitting, four parameter-sine fit 
and seven parameter-sine fit algorithms. The main objective 
is to find a set of parameters to the fitting model (analytical 
expression of a sine signal) that minimizes the sum of the 
squared errors between the model and the sampled data. In 
three parameter sine fit version, the phase, amplitude and DC 
components of an acquired sine wave when the signal rate is 
known and is a multiple linear regression method that requires 
no iterations. In the four parameters algorithm, the rate 
frequency is also calculated which makes the model and 
requires an iterative non-linear least squares procedure in 
order to obtain the best parameters. The convergence of the 
algorithm depends on the initial calculations of the 
parameters. For two channel systems with a common signal 
frequency, as is the case of impedance measurements, all the 
information from both records should be used to obtain a 
better calculate of the common rate of occurrence. This is 
obtained with the seven parameter sine fit which is also an 
iterative algorithm much like four parameter. The main 
difference is that it calculates, in each iteration, the two 
amplitudes, the two phases, the two DC items and the 
common frequency. The effectiveness of this method will be 
demonstrated by smaller Cramer Rao lower Bound (CRLB) 
of the frequency and phase difference when compared with 
the application of two four parameters sine-fit. In the seven 
parameter sine-fit [12], the calculated parameter vector for 
each iteration is 


& = 



T 


(3.8) 

Where A k and B k in equation (3.8), are the in-phase and 
quadrature amplitudes from which the sine amplitudes B k and 
phases 0^ are obtained. Af 1 is the correction that updates the 
calculated common frequency. The iteration process ends 


when the relative rate of occurrence adjustment ~ is below a 


present threshold. After the Iterative part of the algorithm is 
completed and the frequency is determined, the amplitudes 
A k and B k are determined with the three parameter sine fit. 
The calculated parameter vector is obtained from 


Where y is the concatenated sample vector 

y=[Ul.1 Uu ... UlJV 1/2.1 1/2,2 ... Ifcjff 


And 


D li 


,<o _ 


Q[' } 

r, H,) 

0 

0 


qT 


With (co=2Tif') 


qP = 


cos(ct>f ftl ) sin(cot*i) 1 
cos(ct>fj, 2 ) sinfcufj^) 1 

,cos(wf fr . w ) sin(o>f*. N ) 1 


(3.10) 


(3.11) 


(3.12) 


And 


r i‘i_ 

‘k “ 


*t! 


Witll = -1 sin(ffly + Mjfty COStohn). 


IAnJ 


(3.13) 


The calculation of D 3 requires 2N*7 words (in DSP 
execution, each word usually corresponds to a 32-bit long 
single accurate float |D^| D 3 requires an additional 7*7 
words while D Y requires 7 words. Overall, the technique 
requires 17N + 63 words [11] (the 17N part corresponds to 
2N for the samples, 14N for D 3 and IN for FFT used in the 
interpolated discrete Fourier transform to calculate the initial 
frequency. Due to the restricted DSP internal memory and 
the space occupied by the program itself and other internal 
variables, the number of samples N is limited. 


Ey 
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Ey 
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Ey 
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E4.7 

0 

0 

0 
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0 
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(3.14) 

With 
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3.4. Ellipse-fit algorithm 


f,, = Y cos 2 ((ot„).E u = Y cos((/Jt„)sin(wt,,) t 

rial rial 

N N 

£lj = Y C0S(Wt n ).£|,4 = Y COS((/Jf„)«, „. 

n »1 rt =1 

N N N 

£22 = Y sin2 ( Wt n)- E 2J = Y sin ( <ot -) £ 2.4 = Y Sin ( t0t ») a 

rial rta | rial 

A ~ Y *»• 4 = Y (*!•!» + ^2 ") • 

Hal Hal V 

N N N 

£45 = Y cos((or n )a 2n ,£ 4 6 = Y sln(«i n )*2ji A? = Y a2n * 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 


The ellipse-fit algorithm estimates the ellipse parameters that 
best fit the XY pairs of voltages from the two channels. From 
the ellipse parameters, the sine amplitudes, dc items and the 
phase difference are determined [11] [12]. The ellipse-fit 
algorithm is a non-iterative method based on Lagrange 
multipliers. The refined execution, first published requires the 
construction of small matrices (3*3) instead of the N*3 large 
matrix used in their multiplication and inversion and also 
determination of the eigenvectors of a 3 by 3 matrix. Through 
mathematical manipulation the common frequency can be 
removed 


| D r ' D' is a symmetric matrix. Finally for | D ¥ 

JV 

J2 cos(a *t n )u A n 

ra=l 

N 

D sin(ct>t n )Uj „ 

n=l 

N 

Y “l ,n 

11 = 1 


[D (0 ] T y = 


N 

+ « 2 ,n« 2 ,n 

11=1 


JV 


2: 

11=1 




2: 

n=l 


COS(Q>tn )u 2 .„ 

s in (a>t n )u 2 „ 

N 

Y " 2 -r, 

n=l 


(3.20) 


With this method the memory usage is reduced to 3N +63 
(the 3N part corresponds to 2N for the samples and IN for the 
FFT used). For the initial A k and calculates, the 
three-parameter sine-fit is applied to each channel using the 
rate of occurrence obtained with the IpDFT. These two 
algorithms are also executed in the DSP. Samples acquired 
with the proposed system at 1 kHz measurement frequency is 
sampled at 48 kS/s are shown together with the reconstructed 
sine signals using the seven parameter algorithm are shown in 
Fig. 3.6. 


- C,\ 2 , (U2 -C 2 \ 2 ,(u,-C,)(u a -C 2 ) 

\ Dt / \ ^2 / O1D2 

x cos(A^) - sin 2 (A^) = 0, 

(3.21) 


equation (3.21). corresponds to a conic 

F(ui. ui ) = au]+but u 2 +cijj + dih + eu 2 + g = 0, 

(3.22) 


Which describes an ellipse when, 

b 2 - 4 ac < 0. 

(3.23) 

By introducing a scaling constraint q, this condition can be 
set to 4ac - = 1. For conic to correspond to ellipse three 

conditions must be verified Di - 0; D a - 0 and 
cos 2 (DJ < r either a line segment (whenever least one of 
conditions is not true) or an ellipse. From the ellipse 
parameters, the sine parameters are obtained through 

n 1 n __L 

(3 24) 

r 2 «J -be r 2 ae - db 

' ~b 2 -4ac 1 ~ b 2 — 4ac (3 25) 



Fig. 3.6. Acquired samples from two channels with (+) for 
channel 1 and (x) for channel 2 and the corresponding sine-fit 
reconstructed signals (lines). [16] 


And 

m (^)— C3 ~ 6 '> 

_ j a (3,27) 

n, _ Vc 

Therefore the scaling constant q does not need to be 
determined. Also, the sign of parameter a (due of the scaling 
constant, required can be replaced by the sign of parameters a 
scaling parameters a & c can both be negative which 
corresponds also to q < 0). The only parameter that needs 
further calculation is the sign of D x which is unavailable but 
can be retrieved from the rotation direction of the ellipse, if 
the ellipse is constructed clockwise > 0 and D^O if it is 
constructed counterclockwise. The sign of the phase 
difference between two consecutive sample pairs [ 11 ] 
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( 3 . 28 ) 

can be used to determine whether the ellipse is being 
constructed clockwise or counterclockwise. The ellipse center 
is calculated using all the required samples. However, due to 
the presence of noise, some of the consecutive samples may 
give the wrong rotation direction, therefore voting system was 
executed. The sign of the sum of the defined also determines 
the sign of Dj_ [ 11 ] 


sign( A<f>) = -sign 



( 3 . 29 ) 


The ellipse-fit requires 2N + 42 words [11] (2N for the 
samples, four 3 by 3 matrix and two element vectors). With 
the current sine-fit exestuation, the memory needed by the 
ellipse-fit is not a big improvement from the memory needed 
by the seven-parameter sine-fit. The main merit is the total 
number of operation - not only for the matrix construction 
and manipulation but also due to the fact that sine-fit is 
iterative and the ellipse-fit is not, examples acquired at 48 
kS/s with the purposed system at 1 kHz measurement 
frequency are shown together with the reconstructed ellipse 
obtained using the ellipse-fit. Is shown in Fig. 3.7 



Fig. 3.7. Acquired samples from two channels (+) and the 
corresponding reconstructed ellipse (line).[16] 

3.5. GoertzeVs algorithm 

Goertzel algorithm is used to minimize the computation 
cost of DFT by almost a factor of two. It is useful in 
applications that require only a few DFT frequency samples. 

Some applications like frequency shift keying 
demodulation or DTMF, where typically two frequencies are 
used to transmit binary data, the circuit is designed only to 
identify the line for two simultaneous frequencies. Goertzel 
algorithm reduces the complex multiplication for computing 
the DFT relatively by a factor of two to the direct computation 
using the equation (3.30). Goertzel algorithm is derived by 


converting the DFT equation (3.30), into a form of 
convolution which is an equivalent form for the DFT equation 
(3.30). For detailed description the Goertzel algorithm is 
explained by terms of mathematical equations 

The DFT X[k] of a disc ( 3 - 28 ) 

(4.11) 

ce x(n) is given, 

H-l 

sr~* -l Ituik 

X[k] = > x(n)e M k = 0.1.2.. 

..N — 1 (3.30) 

E=fi 


H-l 


= V x(n)V< 

(3.31) 

E = E 


Where W = e^T and W = W sH 


H-l 

v ( 3 - 29 ) 

Zj (4.11) 

n = C 

(3.32) 

H-l 

= ^ x(n)wf _rJ - lk 

(3.34) 

n =0 


= (((Wh* x(0) + x(l> )W N - k + x(2) )W N - k + ■ 

+ x(N — l))WH k (3.35) 


The expression in the equation (3.35) can be expressed as a 
recursive equation, 

y (n) = WH k y (n - 1) -h x(n) (3.36) 

Where y C— IJ = 0, The coefficients will be equal to the 
output of difference equation at time n=N. 

XCfc) = 7 00 (3.37) 

Expressing the difference equation in terms of Z-transform 
and the resulting transfer function can be obtained by 
multiplying the numerator and denominator by 1 — Wjf z; “ 1 , 

Yfe) _J_ 1-WflV 1 _ 1-Wfr" 1 

l(z) = m = 1-W N V 1= 1 - ((Wj+WjJz' 1 - z' J ) = ! _ (2 cos z" 1 - z" J ) ^ 

The realization of Goertzel system is represented below in 
Fig. 3.8 



Fig. 3.8. Realization of Goertzel system 
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We do not compute y (n) for all values of n, but only for n = 
N. It implies that we compute only the recursive part, or just 
the left side in the graphical representation of the realization 
Fig. 3.8. for n = [0, 1, . . . , N], which involves only one real 
and complex product rather than a complex and complex 
product as in a direct DFT , plus one complex multiplication 
to get y (N) = X (k). 

IV. PARTICULAR SAMPLING METHOD 
IMPLEMENTATION. 

The particular sampling method is based on taking signal 
samples in exactly determined time moments, allowing to 
simply calculating the vector of fundamental harmonic of the 
measurement signal. When comparing to DFT, the particular 
sampling method uses only a summation of the collected 
samples, and the obtained two sums determine the orthogonal 
parts of the calculated sinusoidal signal. 

Table 4.1. Example values and eliminated harmonics[4] 


?i 

Eliminated harmonics 

3 

3,9,15,21,27,33,39,45,51,... 

5 

5,15,25,35,45,55,65,75,.. 

7 

7,21,35,49,63,77,.. 

11 

11,33,55,77,99,. 

13 

13,39,65,91,. 

17 

17,51,85,. 

19 

19,57,95,. 

23 

23,69,. 


To determine orthogonal parts of voltage and current 
proportional signals we need to acquire samples and calculate 
two sums for each signal. Each sum items have “+” or “-“signs 
depending on parameter m. For easier implementation, we 
separate items of each sum for those with “+” sign (Rep and 
Inip) and those with “-“ sign (Re^ and Im m ). The main aim 
for this implementation of the particular sampling algorithm is 
determination of the orthogonal parts of the fundamental 
harmonic, but not removing the possibly high number of 
higher harmonics, so the Q was assumed as equal 2 and 
according to Table 4.1, = 3 and q 2 = 5. To correctly 

realize the particular sampling method (to sample at certain 
time moments) it is necessary to assure the sampling 
frequency fs is related to measurement frequency f: 

For assumed values of Q=2, =3 and q 2 =5 the formula 

can be evaluated as below: and this means that for the 
assumed parameters, the sampling frequency should be atleast 
60 times or higher than the measurement frequency ,we need 
to acquire 60 samples during a measurement signal period. 
Using the defined V value, we can express sampling moments 
as a sample number by entering a variable D (given in 
degrees). 


Table 4.2. Sampling schedule for the proposed particular 
sampling method implementation[4] 


Sum 

iteration 

number 

Sampling moment according to (11) and (12) 

Sampling moments expressed as D 
multiple-sample number 

Re component 

Im component 

Rep 

Rcm 

Imp 

IfflM 

1 

-x 

-(x+W) 


1 


15 

2 

x+m 

.f+90+ISO 

30 


45 


3 

-(*+y) 

-(jr+00+y) 


10 


25 

4 

X+180+ — 

a 

I+90+I80+ — 

8 

40 


55 


5 

, 180, 

, 180, 


6 


21 

6 

I+180+ — 

5 

101 

X+90+180+ — 

5 

36 


51 


7 

, ISO ISO, 
~(x+ —+■— 

3 5 

, 180 180, 
-{1490 + y + y) 


16 


31 

8 

ISO 180 
J+1S0+—+ — 

3 5 

r+ 00 + .80 + M 

3 5 

46 


61*1 



The sample numbers given in Table 4.2 above can be 
directly used in software for calculation of Re TJL and Im m 
parts of voltage and current. 


4.1. Particular sampling method implementation test by 
simulation 

The particular sampling method presented will be first 
implemented as a Matlab script and tested. As a reference the 
DFT-based method was used assuming the same sampling 
frequency: all acquired samples containing one period of the 
measurement signal were used for DFT, but selected samples 
were used for the particular sampling method. To create a 
simulation situation similar to the real-life one, the generated 
sinusoidal signal has a white-noise signal (normal 
distribution) and 50 Hz noise added (reflecting Interference 
caused by power lines). 

The Matlab program simulates the experimental setup and 
compares the impedance values that are calculated using the 
DFT algorithm and then simultaneously also using the 
concept of Particular sampling. The program is intended to 
calculate the theoretical phase difference error and the error in 
signal ratio which the relative error in calculating the 
impedance using particular sampling. 

After the simulation the graphs resulting in the phase 
difference error and signal ratio error are obtained. 



Fig. 4.1. Error in magnitude estimation by DFT using matlab 
simulation 
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Fig. 4. 2. Error in phase determination using DFT by matlab 
simulation 
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Noise amplitude 

Fig. 4.5. Noise amplitude versus the error in magnitude by 
DFT and particular sampling using matlab simulation 


Fig. 4.1. and Fig. 4.2. shows the error in determining the 
magnitude and phase respectively using DFT, as a means of 
matlab simulation for 100 iterations 



Number of iterations 


Fig. 4.3. Error in magnitude estimation by particular sampling 
using matlab simulation 



Fig. 4.4. Error in phase determination by particular sampling 
using matlab simulation 

Fig. 4.3. and Fig. 4.4. shows the error in determining the 
magnitude and phase respectively using particular sampling, 
as a means of matlab simulation for 100 iterations. 
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Fig. 4.6. Noise amplitude versus the error in phase difference 
by DFT and particular sampling using matlab simulation 
Fig. 4.5. shows the error in determining the magnitude and 
phase respectively using DFT versus noise amplitude and Fig. 
4.6. shows the error in determining the magnitude and phase 
respectively using particular sampling versus noise amplitude 
using matlab simulation. 

As a result of simulation it can be noted from the Fig. 4.1. 
and Fig. 4.2. that the error in computing the impedance using 
particular sampling is more as compared to the computation 
using DFT. From Fig. 4.3. and Fig. 4.4., the error in 
determining the phase using particular sampling is more in 
magnitude as compared to error by using DFT. For low 
frequencies the particular sampling method is more prone to 
noise as it can be observed from Fig. 4.5. and Fig. 4.6. the 
error rate increase as frequency increase as due to the fact that 
there is some data loss when a part of samples are considered 

4.2. Particular sampling-based impedance measurement 
method implementation in PSoC 

A block diagram and a view of the measurement system 
prototype for experimental evaluation of the particular 
sampling method implementation are shown in Fig. 4.7.a) and 
Fig. 4.7.b). The system consists of 3 parts: a PC computer 
which allows to control the device and visualize [4]. 

Thanks to the use of PSoC, the number of items was 
reduced to the minimum. The used PSoC generation 
represents microcontrollers with relatively low processing 
power and small RAM memory. The prototype was built 
using a CY8C29566 chip with 2 kB SRAM memory. The 
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sinusoidal excitation signal applied to the calculated 
impedance is produced using the DDS method with the aid of 
a D/A converter DAC1, on the basis of sine samples placed in 
PSoC’s RAM memory (GENbuf). TIMER1 creates a clock 
signal which controls generation and acquisition using the 
microcontroller clock. The DAC1 output signal before 
application to the impedance under test is first filtered in a 
low-pass filter removing unwanted stair-steps. 

a) b) 



Fig. 4.7. Block diagram (a) and a view of prototype (b) of the 
particular sampling method PSoC implementation. [4] 


4.3. Implementation in form of Labview application using 
DAQ card 

In this method, a DAQ card is used to produce the 
excitation signal that is the voltage signal that is used to find 
the impedance .The FFT of the signal is taken and displayed 
as phase and amplitude spectrum. The excitation signal is 
passed through the circuit under test and the current through 
the circuit is observed and voltage across it is calculated. The 
resultant values are then transformed to FFT and its phase and 
amplitude spectrum is observed, thus the impedance can be 
observed by dividing the amplitudes and subtracting the 
phases. The particular sampling method is based on taking 
signal samples in precisely determined time moments, 
allowing to simply calculate the vector of fundamental 
harmonic of the measurement signal. When comparing to 
DFT, the particular sampling method uses only a summation 
of the collected samples, and the obtained two sums 
determine the orthogonal parts of the measured sinusoidal 
signal. 

The front panel of the signal generator developed by the 
labview application is shown below in Fig. 4.8. 



Fig. 4.8. Front Panel of labview application for signal 
generation 


The application provides functionality to the user to 
generate the different kinds of waveform, change the 
amplitude and frequency of the signal. The Cycles per buffer 
is the number of cycles viewed on the display and most 
preferably it is set to 1. 



Fig. 4.9. Front panel of labview application for sample 
acquisition 


The application provides the functionality to change the 
sampling frequency to obtain different sampling ratios which 
are useful for experiments and is shown in Fig. 4.9. There are 
two channels, whose signals are acquired using the DAQ as 
shown in Fig. 4.10. to analyze the error in the method devised 



Fig. 4.10. DAQ card NI USB - 6251 used for experiment [19] 
Experimental setup 

An excitation signal is generated using the DAQ NI USB 
6251, to the circuit as shown in Fig. 4.11. 



Fig. 4.11. Circuit diagram showing experimental setup 
(DIPTRACE software) 
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This is actually a bridge setup where the impedance to be 
found is connected to the unknown terminal and the known 
impedances are connected to other terminals. The excitation 
signal is sine wave and the signal are back fed to DAQ and 
processed by the lab view application to compare the phase 
difference and the ratio of the signal amplitude using FFT and 
particular sampling. Since the application provides the ability 
to the user to change amplitude, frequency, cycles per buffer. 
The experiment procedure is as follow 


•Varying the amplitude of the signal for 3 different 
volts as 1 volt, 3volts and 5 volts. 

•Varying the sampling ratio to 60, 80, 100- by 
changing the sampling frequency 
•Adding additional component in the impedance to 
be calculated 

•Incrementing the sample index that are used to 
obtain the particular samples 
•Varying the frequency of the signal 
Experiment 1- Using only resistor in all positions 


Resistors of value 10K are connected in all terminals. The 
amplitude of the signal is varied from 1, 3 and 5 Volts and the 
frequency varied from 0.1, 1, 10, 100, 500, 1000, 2500, 5000 
Hz and also the sampling ratio that is the ratio of sampling 
frequency to actual frequency is varied from 60 80 and 100. 
All of which has been carefully depicted in graphical 
representation. 



-at sampling ratio 
60 

-at sampling ratio 
80 

-at sampln ratio 
100 


1000 2000 3000 4000 SOOO 6000 
Frequency 


a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 



—at sampling 
ratio GO 

—at sampling 
ratio 80 

— at sampling 
ratio 100 


F r equency 



—at sampling 
ratio 60 

—at sampling 
rate 80 
at sampling 
ratio 100 


Frequency 



Frequency 



Frequency 


C) 



0 1000 2000 3000 4000 5000 6000 

Frequency 

Fig. 4.12. Error in phase difference 


Fig. 4.13. Error in signal ratio 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts. 


The experimental results for the error in phase difference and 
signal ratio for different voltages 1 V, 3 V and 5V, and 
different sampling ratios 60, 80 and 100 are shown in Fig. 
4.12 and Fig. 4.13. 

Experiment with an additional capacitor connected parallel 
to resistor The resistor of value 10K and a capacitor of 100 
pF are connected parallel across the terminal for impedance 
measurement 

a) 



Frequency 
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b) 

072 



Frequency 

C) 

0 5 t- 
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Frequency 



at sampling rat© 
60 

—at sampling rat© 
80 

— at sampling rat© 
100 


Fig. 4.14. Error in phase difference 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 


C) 

0.09 



Frequency 


Fig. 4.15. Error in signal ratio 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts. 


The experimental results for the error in phase difference 
and signal ratio for different voltages 1 V, 3 V and 5V, and 
different sampling ratios 60, 80 and 100 are shown in Fig. 
4.14. and Fig. 4.15. 

Experiment 3 - With additional two capacitors connected 
parallel to resistor. 

The resistor of value 10K and two capacitors of values 100 
pF and 100 nF are connected parallel across the terminal for 
impedance measurement. 



Frequency 


b) 


b) 



Frequency 


045 



0 1000 2000 3000 4000 5000 6000 
Frequency 


at sampling 
ratio 60 

at sampling 
ratio 80 

at sampling 
ratio 100 
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The experimental results for the error in phase difference and 
signal ratio for different voltages 1 V, 3 V and 5V, and 
different sampling ratios 60, 80 and 100 are shown in Fig. 
-at samping 4.16. and Fig. 4.17. 

ratio 60 Experiment 4 - The sample index are incremented by 1 

a) 

-at sampling 
ratio 80 


1000 2000 3000 4000 5000 6000 
Frequency 


-at sampling 
ratio 100 


Fig. 4.16. Error in phase difference 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 


a) 



at sampling 
ratio 60 

at sampling 
ratio 80 


ratio 100 


1000 2000 3000 4000 5000 6000 
Frequency 


b) 



-at sampling 
ratio 60 

-at sampling 
ratio 80 

•at sampling 
ratio 100 


Frequency 


C) 
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at sampling 
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at sampling 
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Frequency 

Fig. 4.17. Error in signal ratio 
Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 
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Frequency 
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Fig. 4.18. Error in phase difference 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 
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b) 


c) 


0.2 
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a) 


Frequency 

Fig. 4.19. Error in signal ratio 
Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 


The experimental results for the error in phase difference and 
signal ratio for different voltages 1 V, 3 V and 5V, and 
different sampling ratios 60, 80 and 100 are shown in Fig. 
4.18. and Fig. 4.19. 

Experiment 5 - The sample index incremented by 3 

a) 
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Fig.4.20. Error in phase difference 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 
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Fig.4.21. Error in signal ratio 
a) Using 1 Volt, b) Using 3 Volts, c) Using 5 Volts 
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The experimental results for the error in phase difference and 
signal ratio for different voltages 1 V, 3 V and 5V, and 
different sampling ratios 60, 80 and 100 are shown in Fig. 
4.20. and Fig. 4.21 

V. RESULTS 

It can be seen that, there are many advantages and 
disadvantages using particular sampling method over the 
most accepted DFT algorithms. As from the the experimental 
results of the particular sampling using the Lab-view 
application using DAQ card shows the main inference as that 
the error in computing the impedance is affected by the signal 
frequency and the voltage doesn’t have any vast influence in 
the computation of the impedance 

The mean value of error in computation of phase difference 
is plotted as a function of frequency is shown in Fig. 5.1. and 
it shows, increase in the frequency resulted in the error rate up 
to a certain extent 



Frequency 

Fig. 5.1. Graphical result of value of error in phase difference 



0 1000 2000 3000 4000 5000 6000 

Frequency 

Fig. 5.2 . Result of mean value of error in signal ratio 

The Fig. 5.2. shows the variation of the signal ratio error as 
a result of frequency change and it is evident from the graph 
that the errors increase with the increase in frequency. 

The Fig. 5.3. shows the average percentage of error in 
computation of impedance using particular sampling method 
and it can be seen that the overall error in percentage varies 
from 2% in low frequencies to an average of 6% in high 
frequency. 



Frequency 

Fig. 5.3. Graphical result percentage of error for various 
frequencies 

Advantages of particular sampling method 

For N samples M 2 complex multiplications and N(N-l) 
complex additions are needed for DFT algorithm So in our 
experiment (particular sampling takes minimum 60 samples). 
Total number of calculations using DFT will be 7140 
complex multiplications and additions. For Particular 
sampling we need only additions and the total additions will 
be 4N. So total calculations using particular sampling is 240 

So when the time saving is considered, with the 
implementation of the particular sampling algorithm we will 
be saving a large amount of time and memory as the same for 
the same result as compared to DFT algorithm with a certain 
amount of error. As precisely there is saving of 96.3% time 
and memory if employ the particular sampling for the same 
purpose using DFT. As per the experimental results the 
maximum error in the calculation of impedance measurement 
is 6%. 

It should be noted that the saving of time can be achieved in 
the loss of accuracy of the measurement with an average 
standard deviation of 0.041 

VI. CONCLUSION 

The use of the particular sampling method for impedance 
measurement allows to meaningfully decrease the number of 
samples and the calculation effort. For example, for the 
proposed realization only 8 samples for each part (Re, Im) of 
each signal (voltage, current) are acquired and then added. On 
the other side, the accuracy of the method as well as the 
immunity to noise is lower than in the DFT-based impedance 
measurement. The obtained accuracy is acceptable in case of 
measurement in the field and the resulting simplification of 
the device can open a new application area for impedance 
measurement (e.g. smart impedance sensors). The developed 
system is a prototype, which has made possible experimental 
verification of the method using particular sampling. The 
system has created a base for development of a simple, small 
impedance measuring application. 

From the experimental results and comparing the DFT and 
particular sampling algorithms, it can be noted that DFT is 
more reliable and is accurate and more time consuming than 
the other methods. The memory needs to process a lot of 
mathematical data is more so it is the most possible 
disadvantage of DFT. While on the implementation of 


31 


www.erpublication.org 


































Analysis of Impedance Measurement Implementation Using Particular Sampling (By Lab View and Matlab) 


particular sampling method we have a large saving of time 
and memory more that 90% and only a minute loss of 
accuracy of 6%. So the time is saved on expenditure of 
accuracy. The particular sampling methods provide an easy 
and not tiring method for impedance measurement and can be 
used for various domains of consumer application especially 
commercial. 
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Fig. 3.6. Acquired samples from two channels with (+) for channel 1 and (x) 

for channel 2 and the corresponding sine-fit reconstructed signals (lines). 

Fig. 3.7. Acquired samples from two channels (+) and the corresponding 
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Fig. 3.8. Realization of Goertzel system 
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simulation 
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particular sampling using matlab simulation 
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Fig. 4.7. Block diagram and a view of prototype 
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Fig. 4.9. Front panel of labview application for sample acquisition 
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Fig. 4.11. Circuit diagram showing experimental setup (DIPTRACE 

software) 
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Fig. 4.14. Error in phase difference (experiment 2) 
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Fig. 4.16. Error in phase difference (experiment 3) 

Fig. 4.17. Error in signal ratio (experiment 3) 

Fig. 4.18. Error in phase difference (experiment 4) 

Fig. 4.19. Error in signal ratio (experiment 4) 

Fig. 4.20. Error in phase difference (experiment 5) 

Fig. 4.21. Error in signal ratio (experiment 5) 

Fig. 5.1.Graphical result of value of error in phase difference 

Fig. 5.2. Result of mean value of error in signal ratio 

Fig..5.3. Graphical result percentage of error for various frequency 
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